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The stability of current atomic clocks is limited by the quantum noise of the atoms. To reduce 
this noise it has been suggested to use entangled atomic ensembles with reduced atomic noise. 
Potentially this can push the stability all the way to the limit allowed by the Heisenberg uncertainty 
relation, which is denoted as the Heisenberg limit. In practice, however, entangled states are often 
more prone to decoherence, which may prevent reaching this performance. Here we present an 
adaptive measurement protocol that in the presence of a realistic source of decoherence enables us 
to get near Heisenberg limited stability of atomic clocks using entangled atoms. The protocol may 
realize the full potential of entanglement for quantum metrology despite the detrimental influence 
of decoherence. 
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Some of the most precise measurements in physics are 
made with atomic clocks. Very accurate time measure- 
ments are used in anything from GPS-systems to gravi- 
tational wave detectors. Currently one of the main lim- 
itations to the stability of atomic clocks comes from the 
quantum noise of the atoms, which leads to the stan- 
dard quantum limit (SQL) where the resolution scales as 
I/Vn with N being the number of atoms [l][2]. To over- 
come this noise it has been suggested to use entangled 
states with reduced atomic noise [SHZ]- Ultimately this 
may lead to a stability at the Heisenberg limit where the 
resolution scales as 1/N, and recently the first proof of 
principle experiments have demonstrated these concepts 
experimentally [8Hl3]. In practice, however, entangled 
states are often more prone to decoherence, and to fully 
asses the advantage it is essential to study the perfor- 
mance in the presence of decoherence [H]. In Ref. [5 
it was proven that entanglement can be used to improve 
the long-term stability of atomic clocks in the presence 
of the dominant practical source of decoherence, but the 
improvement identified was rather limited. Here we show 
that it is possible to obtain a large improvement in the 
resolution of the atomic clock by combining entangle- 
ment with an adaptive measurement protocol (inspired 
by Ref. [E]). With our adaptive measurement protocol 
the entangled states are not more sensitive to the deco- 
herence than disentangled states. As a consequence the 
long term stability of the atomic clock can be improved 
almost to the Heisenberg limit even in the presence of 
decoherence. 

An atomic clock operates by locking a local oscilla- 
tor (LO) to an atomic transition via a feedback loop. 
The feedback is based on a measurement of the LO fre- 
quency offset Sluq compared to the atomic transition. The 
frequency offset is typically obtained through Ramsey 
spectroscopy [16] . Ramsey spectroscopy consists of three 
parts; first the atoms are illuminated by a near-resonant 
7r/2-pulse from the LO followed by the Ramsey time T 



of free evolution, and finally another near-resonant 7r/2- 
pulse is applied to the atoms. During the free evolu- 
tion the LO acquires a phase S(j)o — SuioT relative to the 
atoms. Due to the 7r/2-pulses this phase can be mea- 
sured as a population difference between the two clock 
levels. SujQ can thus be estimated from the measurement 
and used for a feedback that steers the frequency of the 
LO to the clock-frequency. The resolution of the clock 
will improve with T since a longer T improves the rel- 
ative sensitivity of the phase measurement. For current 
atomic fountain clocks, T is limited by gravity and can 
hardly be varied. Here on the other hand we consider 
trapped particles, where T can be increased until it is 
limited by the decoherence in the system. The long term 
stability thus depends on the nature of the decoherence. 
To take this into account Ref. [14^ considered dephasing 
of single atoms to be the dominant noise in the system. 
For this model Ref. [T3] showed that entanglement can 
not improve the stability of atomic clocks considerable 
(although an improvement is possible for non-Markovian 
noise [U]). A more realistic model of the decoherence 
was described in Ref. [S] where the primary noise source 
is the frequency fluctuations of the LO [TH] . In this work 
a small improvement in the long term stability, scaling as 
^ iV^/^, was identified for entangled atoms. Here we use 
the same realistic model of the LO noise and show that 
entanglement and adaptive measurements may improve 
the performance and give near Heisenberg limited atomic 
clocks. 



We consider an ensemble of N two-level atoms, which 
we model as a collection of spin-1/2 particles. The to- 
tal angular momentum is J = X)i=i ^i where Si is the 
spin-vector of the i'th atom. The angular momentum 
operators Jx,y,z gives the projections of J on the x, y 
and z-axis. The atoms are initially pumped to have a 
mean spin along the z-axis, (Jx) — (Jy) — 0. After 
the Ramsey sequence the Heisenberg evolution of Jx , Jy 
and Jz is Ji = Jx,J2 = sin(J(/)o)Jy — cos((S0o)^z and 




FIG. 1. (Colour online) The atomic state just before the 
measurement of J^ for (A) uncorrelated atoms, (B) moder- 
ate squeezed atoms and (C) highly squeezed atoms. Highly 
squeezed atoms introduce significant extra noise in the mea- 
surement through a noise term oc 5(j!)AJz/|{Jz)| in the phase 
estimate. The adaptive protocol seeks to rotate the atomic 
state to have mean spin almost along the J^-axis thus reduc- 
ing this noise. 



J3 = cos{6(j)o)Jy + sm{S<j)Q)Jz. At the end of the Ramsey 
sequence J3 is measured and used to estimate d(j)o. The 
Jy term in J3 results in the so called projection noise 
in the phase estimate ~ AJj,/|(J2)|. For uncorrelated 
atoms AJyAJr^ — {Jz)/'^ ~ N/A and the projection noise 
causes the stability of the clock to scale as ^ 1/\/N. 
For a spin squeezed state |19j the variance of Jy is re- 
duced to obtain a better phase estimate. Such a spin 
squeezed state is depicted in Fig. [T] which shows how 
the noise of the spin squeezed state looks like a "flat ba- 
nana" lying on the Bloch sphere. The more we squeeze, 
the longer and more narrow the banana is and significant 
extra noise is added to the mean spin direction. For a 
phase measurement based on a direct measurement of J3 
this gives an additional noise term ^ 5(I)qAJz/\{Jz)\ in 
the phase estimate. This extra noise limited the perfor- 
mance in Ref. [5] if strongly squeezed states were used 
such that only a limited improvement could be identi- 
fied. We avoid this problem by using an adaptive scheme 
with weak measurements to make a rough estimate of 
(500 and then rotate the spins of the atoms such that 
the mean spin direction is almost along the y-axis. The 
fiat banana depicted in Fig. [T] will afterwards lie in the 
xy-plane and this will eliminate any noise from AJ^ in 
subsequent measurements. Having eliminated the noise 
from AJz we can allow strong squeezing in AJy and ob- 
tain near Heisenberg limited stability. 

To be concrete our weak measurements is based on the 
strategy developed and demonstrated in Ref [101 IHl EO] 
where a light field dispersively interact with the spin 
and is subsequently measured. This is described by a 
Hamiltonian Hint = —XiJs-^i where xi is the interac- 



tions strength and Xi is the canonical position operator 
of the light PTH23] . The measurement results in a rota- 
tion around J3 described by the rotation matrix R3(ni), 
where 11 = fliXi. ili = XiA*i is the measurement 
strength and fj,i is the measurement time. The canon- 
ical momentum operators of the light before A and after 
P, the interaction are then related by Pi = Pi — ili J3. 
Pi is measured using homodyne detection [24] and the 

phase is estimated as Scjyf — ~ '"- \ where the factor /3i 

is found from minimizing (((5(/)o — S(j)f)'^). Based on the 
phase estimate we rotate the spin of the atoms around 
Ji in order to compensate for the extra noise added by 
the spin squeezing. This is described by a rotation ma- 
trix Ri(J0f ). The process can be iterated such that after 
n—1 weak measurements the Heisenberg evolution of the 
original operators ( Ji, J2, J3) is: 

i I = Ri(50^_i)R3(n„-i)...Ri WDR3(ni)( i2 1 (1) 
kJsJ „ Vis/ 

The final measurement is assumed to be a projective 
measurement and the final phase estimate Scj)^ is thus 






The factors of (3i in the phase estimates 

are found be minimizing ((i50o~X]i=i ^4'Tf') with respect 
to Pi after each measurement. 

We will now show semi-analytically that the measure- 
ment strategy in equation Q allows for near Heisenberg 
limited stability. For simplicity we set all Pi — 1. After 
a series of n adaptive measurements the estimate 54>'^ of 
the initial phase of the atoms (^0o) is 5(t>'^ — X]r=i ^^t- 
6(j)'^ is then used for the feedback on the LO. We as- 
sume that negligible time is spent on everything but 
the free evolution period T. After a time tk = kT the 
frequency correction Auj{tk) = —a6(t>e{tk)/T is applied 
to the LO where a sets the strength of the feedback 
loop. The frequency offset of the LO at time tk is then 
Su}{tk) = Sujo{tk)+J2i=i ^'^i'ti)i where Sujo{tk) is the fre- 
quency fluctuation of the unlocked LO and Acu {ti) is the 
frequency correction applied at time ij. The mean fre- 
quency offset after running for a period t = IT {I ^ I) 
is 



fc(r) = ] y: 



k=l 



S4>{tk) - S(j)e{tk) 

T 



(2) 



where S(j){tk) is the accumulated phase of the LO during 
the free evolution between tk-i and tk and S(j)e(tk) is the 
estimate of this phase (see supplementary information) . 
The measure of the long term stability of the atomic clock 
is: 



a,{r) = {iSuir)/Lo)y^' 




T 



(3) 

(4) 



We initially assume that the phase offset of the unlocked 
LO (500 is due to phase fluctuations in the LO with a 
white noise spectrum. Later we will also consider the 
case where the fluctuations have a 1// spectrum. For 
white noise we have ((5(/)o) = ^T {{6cf>o) = 0) where 7 is a 
parameter characterizing the fluctuations of the LO. In 
the limit of a ^ 1, the phases are uncorrelated such that 
cr^(T) = ^7/™2(((50o-,S(/,'=)2)/^T)i/2. This expression 
shows that for fixed 7, r the stability of the clock only 
depends on how precisely we can estimate (5(^o- 

We now turn to the adaptive measurement strategy to 
estimate d4>o- After j weak measurements the difference 
between the true phase and the estimated phase 6^j is: 



6^, 
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HI (5) 



Using equation M to get an expression for Ji/)^ and the 

fact that ^$j-i — 5<f>o — X]i=i Hi^ '^^ can express the 
phase error as 

where we have assumed (5$j_i<Cl. The first term in ([6]) 
gives a contribution ^ (5$j_iAJj,/( J^) to <Tj{t) from the 
noise in the mean spin direction as previously mentioned. 
Note that this term is proportional to the phase estima- 
tion error at the previous measurement stage, since it 
depends on how well the banana in Fig. [l] is rotated into 
the xy-plane. For a useful adaptive protocol (5$j_i gets 
smaller for growing j and the noise that enters through 
AJz is reduced. The last terms in ^ gives the noise 
from AJj,, the accumulated back action of the previous 
measurements {{SJ^j}) and the noise from the incoming 
light in the measurement (AP? = (Pf))- 

The stronger a measurement is, the less noise is added 
through AP'^/n'^ since the measurement is more pre- 
cise. Any imprecision APf^JVlf^, from previous mea- 
surements is contained in 5$j_i and is corrected for in 
the subsequent stages of the protocol, which essentially 
estimate how well we corrected the phase in previous 
measurements. This means that we can initially work 
with weak measurements, which only give a rough esti- 
mate since later stronger measurements correct for the 
imprecision in the initial measurements. 

The accumulated back action noise SJ^^j originates 
from the disturbance caused be the measurements. The 
measurements adds noise in Ji , J2 , which is mixed into 
J3 when the atomic state is rotated to have mean spin al- 
most along the y-axis. From equation (IT]) the dominant 
terms in SJ^j is found to be 5Jjj ^ X]i=i Ht^i^iJx 
(see supplementary information). It is seen that the 
stronger a measurement is, the more noise is added to the 
stability. However for a useful adaptive protocol S(j)'i gets 
smaller for growing i, which again means that the i'th 



measurement can be stronger than the previous (z— l)'th 
measurement without adding more noise to the stability. 

Above we have argued that we can suppress the noise 
terms originating from AJ2, APj, and SJ^j using an 
adaptive protocol with weak initial measurements. A 
remaining question is how well this suppression work. 
This is considered in detail in the supplementary infor- 
mation. To be specific, we consider spin squeezed states 
ofthcformlV'(K)) =7V(K)X;„(-l)"e-(™/'')'|TO), where 
\m) are eigenstates of Jy with eigenvalue m, A/'(k) is a 
normalization constant and the sum is from —J to J 
where J = N/2 is the total angular momentum quantum 
number, k gives the width of the state and thus deter- 
mines the degree of squeezing. Uncorrelated atoms are 
well approximated by k = iV^'^ while highly squeezed 
states are obtained for k -^ 1. We furthermore as- 
sume that the probe light has vacuum statistics such that 
(P) = (X) = (XP) = and (X^) = (P^) = 1/2. As an 
upper limit of the stability wc find that for n ^ 3 \og{N) 
weak measurements, using a spin squeezed state with 
K ~ log \/N + 2 and choosing a measurement strategy 
with D,i ^ ]\[-^+'^/i"-+^) we can suppress other noise terms 
so that the measurement is eventually limited by AJj,. 
ct^(t) will then be -- (2/iV + \ogVN/N)/y/Yr (in units 
of ^-iJTuP-) for TV > 1. It is seen that for N = 10*^ the 
upper limit of the stability will differ from the Heisenberg 
limit, which is <J-^{t) = [1 / N) / \f-{T by a factor of ~ 5. 

The above upper limit shows that we can be near the 
Heisenberg limit, but to get a better estimate we have 
made a numerical minimization of a-^ in order to find the 
optimal resolution. Wc have simulated an atomic clock 
with a LO subject to white Gaussian noise for atom num- 
bers in the range N = 100 to A^ = 10^ (for details see 
supplementary information). For N < 1000 we simu- 
late the full quantum evolution (we denote this as 'full 
quantum simulation') and for N > 1000 we approximate 
the probability distributions of Jx,y,z with Gaussian dis- 
tributions with moments calculated from |'0(k)) in the 
approximation A^ ^ 1 (we denote this as 'Gaussian sim- 
ulation'). 

An example of the results are shown in Fig. [2] which 
demonstrates that there is a significant improvement by 
using a spin squeezed state compared to uncorrelated 
atoms. Also with the adaptive protocol, the Ramsay 
time can be as large for highly entangled states as for 
disentangled state and there is thus no difference in the 
relevant coherence time. Furthermore the adaptive pro- 
tocol allows longer interrogation times 7r < 0.3 than 
the conventional protocol 7T < 0.1. This is because the 
adaptive protocol can determine phases < tt while the 
conventional protocol begin to give ambiguous results for 
phases ~ 7r/2. 

We have numerically minimized cr^(T) in both the de- 
gree of squeezing, the number of weak measurements, 
the Ramsey time, and the strengths of the measure- 




FIG. 2. (Colour online) Stability as a function of the Ramsey 
time (77) for A^ = 10^. D, V are the best non-linear protocol 
of Ref. [5] while o,a are the adaptive protocol discussed in 
this article. The adaptive protocol allows for 'yT ~ 0.3 while 
the conventional protocol of Ref. fSl only allows for 'yT ~ 
0.1. V, A correspond to uncorrelated atoms while o, D are the 
ideal choices of squeezing in the adaptive and conventional 
protocols respectively. For the adaptive protocol the figure 
shows how the probabilistic nature of the errors that limits 7T 
makes the stability jump back and forth for 'yT > 0.3 before 
it is definitely diminished (see supplementary information). 




FIG. 3. (Colour online) Optimized stability of an atomic 
clock for a LO subject to (a) white noise and (b) l//-noise. 
o, D, A, V are the results for the full quantum simulation while 
• ,B,A,T are the results for the Gaussian simulation. The 
Gaussian simulation has been extended down to TV = 100 
and was found to give more or less identical results to the full 
quantum simulation. o,» (A, A) are the adaptive scheme and 
D, ■ (V,T) are for the conventional protocol with (without) 
entanglement. The dotted lines are the analytical results and 
the solid line is the Heisenberg limit for the maximal Ramsay 
time 7T = 0.3 (a) and -fT = 0.2 (b). 



ments. Fig. IsFa) shows the result of the optimization for 
both the adaptive protocol and the conventional protocol 
with/without squeezing. The adaptive protocol gives a 
significant improvement compared to using uncorrelated 
atoms resulting in near Heisenberg limited stability. The 
numerical calculations also agree nicely with the analyt- 
ical calculations (see supplementary information for de- 
tails). As noted above the adaptive protocol allow for 
a longer Ramsey time than the conventional protocol, 
which gives an improvement of roughly a factor 1.6 for 
uncorrelated atoms (see Fig. [2]) . Using the adaptive mea- 



surement protocol thus leads to an improvement even 
without entanglement. 

So far we have assumed white noise in the LO. In prac- 
tice the noise of the LO is however more likely to have a 
nontrivial spectrum like l//-noise. We have therefore re- 
peated the numerical simulations with the LO being sub- 
ject to l//-noise for which {uj{f)iu{f')) = S{f + fh^/J. 
The details of the locking of the LO to the atoms through 
a feedback is described in the supplementary informa- 
tion. We have performed the numerical optimization in 
the same parameters as for white noise. The results are 
shown in Fig. p[b) and shows that the improvement ob- 
tained using the adaptive scheme with correlated atoms 
persists also for 1// noise. Again near Heisenberg lim- 
ited stability is obtained using the adaptive protocol. 
The longer Ramsey time of the adaptive scheme com- 
pared to projective measurements gives an improvement 
of roughly a factor 1.3 for uncorrelated atoms. 

In conclusion we have developed an adaptive measure- 
ment protocol which allows operating atomic clocks near 
the Heisenberg limit using entangled spin squeezed en- 
sembles of atoms. These results clearly demonstrates 
that entanglement can be an important resource for 
quantum metrology. Importantly our results are ob- 
tained under realistic assumptions where we account for 
the dominant source of noise in practice. We find that 
in this situation wc can gain nearly the full potential 
of entanglement estimated without accounting for deco- 
herence. Furthermore the adaptive protocol allows for a 
higher Ramsey time, which gives an improvement even 
for uncorrelated atoms. 
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SUPPLEMENTARY INFORMATION 

This supplementary material to our article " Near Heisenbcrg limited atomic clocks in the presence of decoherence" 
consists of two parts. In the first part we will go through the adaptive protocol in detail and describe the various 
noise terms that effect the stability of the clock. We show semi- analytically that the adaptive protocol gives near 
Heisenbcrg limited clocks with a stability scaling at most like log N/N. The second part describes the details of our 
numerical simulations of an atomic clock. We describe the subtleties of simulating an atomic clock running for a long 
but finite time and show how the feedback on the LO effectively locks it to the atoms. 

SCALING OF THE STABILITY 

In the article we describe how the adaptive protocol consists of a series of weak measurements of the atomic spin 

with intermediate feedback on the atoms. After a weak measurement we estimate the phase of the LO relative to the 

ftp ^ f 

atoms as b(^\ = ' - ' where P^ is the canonical momentum operator of the light after the interaction with the atoms 

and Vli is the measurement strength. To simplify the equations we set Pi = 1 throughout this section. Note that this 
choice is not ideal and the true performance of the clock will thus be better than what we estimate here. In order to 
go through the details of the adaptive protocol we consider a Ramsey sequence where the LO acquires a phase (5(/)o 
relative to the atoms during the free evolution. 54)q is subsequently estimated through n—1 weak measurements with 
intermediate feedback on the atoms and a final projective measurement. We define the operators 

■k = 4, J2{0) = sm{e) Jy - cos{9)J„ .h{e) = sin(0)J, + cos{e)Jy (7) 

such that the outcome of the Heisenbcrg evolution of J^, Jy and Jz during the Ramsey sequence is Ji, J2(<50o) £^nd 
J-^{5(J)q). In the Heisenbcrg picture the first weak measurement results in a rotation of the atomic spin described by 
a rotation matrix R3(Hi), where H = JliXi as described in the article [Xi is the canonical position operator of the 
probe light). The subsequent feedback on the atoms is a rotation described by the rotation matrix Ri(50^) where 
5(j)\ is the phase estimate based on the measurement result. Thus the Heisenbcrg evolution of the operators after the 
first weak measurement and subsequent feedback is: 

Ri(<50?)R3(ni)(j2Wo) (8) 

cos Hi —sin Hi 

cos(50^sin(IIi) cos((5(/)f)cos(Hi) -sin((5</-f) | ( J2{5<i)o) | (9) 

Vsin((50!j;) sin(IIi) s\n{5(f>'i) cos{t\i 

Ji \ f^Jia 






J2((50o-'50f)| +(<5J2,2| (10) 

where the operators 

5.h,2 = (cos(ni)-l)Ji-sin(ni)J2((50o) (11) 

5J2,2 = cos((50^)sin(IIi)Ji+cos(,5(/.^)(cos(ni)-l)J2(<5(/'o) (12) 

ska = sin((5(/.^)sin(ni)Ji+sin((50^)(cos(ni)-l)J2((50o) (13) 

describe the noise due to the back action of the measurement. Note that the back action also effects J3 even though 
J3 is conserved and hence unaffected during the measurement. This is because the rotation during the feedback mixes 
back action noise into J3. 

The process is now iterated such that the Heisenbcrg evolution after n—1 weak measurements and subsequent 
rotations is 

lJl\ I J^ \ /SJl.n\ 

(14) 



•JA 


/ 


Ji 


'&- 


= h {S4>Q - 




jJ. 


\MS(t)o 






where the iterative expressions of 5Ji,j-i,SJ2.j-i,SJ3j-i are 
SJij = (cos(nj_i)-l)Ji-sin(nj_i)J2 




cosfn 



J-1 



SJ2,j — COS 



SJ: 



3,3 



■J_i) sin(n,_i) Ji+ cos(<5(/.^_i)(cos(n,_i)-l) J2 | 
+ cos(5(?!)^_J cos(nj_i)5J2j_i-sin((50^_J(5J3j_i 
= sin((50j_i) sin(nj_i) Ji+ sin((50J_i)(cos(flj_i)-l) J2 | 



)'5Ji,,-i 


-~sm{n^_i)SJ2,j-i 


(15) 




+ cos(50^_i) sin(nj_i)(5 Ji j-_i 


(16) 




+ sm{S(t)'j_i) sin(nj-i)5 Ji j_i 





sm 



tij_i)cOs{tlj^l)SJ2^j-l+COs{S4}j_i)dJ3j^ 



(17) 



with 5Ji.i = SJ2,i — SJ^i = 0. In the final projective measurement we measure J^^n and obtain a phase estimate 
d(j)^ ~ Js^n/iJz)- Our final estimate of (50o is then X]i=i ^4>i and the difference between this and the true phase is 

-<5C. (18) 

Using equation (14) we can express this using the previous phase estimation error ((5$„_i) and the last measurement 



<5$„ = 


= (500 - 


1=1 


= Hq - 


ri-l 



^$. 



^<5$„_i(l- JJ(J,))- Jy/(J,) - W3,„/(J.) 



(19) 
(20) 



where (5$,i_i = (50o ^ X]i=i '^'^i ^^^ we have assumed (5<i>„„i <Cl to expand the sine and cosine. Below we will argue 
that for white noise and in the limit of a weak feedback {a <^ 1) we can determine the stability of the LO from looking 
at the phase error ((5$^)^/^ for each Ramsay sequence independently. We shall therefore now explore the limitations 
arising from the noise term in this expression. 



Equation ( 20 1 shows three kinds of noises that limit the stability. The first term is due to the uncertainty in J^ 



(A J^) and is proportional to (5<i>„_i. This is the noise that limits the performance of the conventional protocol in Ref. 
[5] and the idea of the adaptive protocol is to eliminate this noise by minimizing (5$„-i through many measurements 
and feedback. The second term is the noise in Jy^ which is the term we want to decrease by squeezing the atomic 
spin. The final term is the accumulated back action of the measurements. The weaker the measurements are the 
smaller the contribution from this term will be. 

We now analyze the various terms in ((5<&^i) in more detail. The estimated phase in the j'th weak measurement is 



J0j' = (sin(J$j_i)J^ + cos((5$j_i)Jj, + SJ^.j - Pj/f^j) /(J.). 



(21) 



Using equation (15)-(17l and equation (20l-(21) we can express the dominant contributions to the stability due to 
the noise in J, as 



(<5<^g)((l-J,/(J,))2") 
+2(<50o(sin(J0o) - '5(/)o))((l-i./(i.))'""'i./(iz)) 
+ ((sin(5</)o) - 5<t>of){{l-JrJ{l)f''-^Jll{.h?) 



(22) 
(23) 
(24) 



Here "dominant" refers to decreasing slowest with N . The above expressions are independent of the atomic state but 
we will now focus on a specific type of states in order to treat the system in more detail. 

We consider spin squeezed states of the form |?A(k)) — A/'(k) X]m(~l)'"^~*''"^''^ 1"^)' where \m) are eigenstates of 
Jy with eigenvalue to, M{k) is a normalization constant and the sum is from —J to J where J — N/2 is the total 
angular momentum quantum number, n gives the width of the state and thus determines the degree of squeezing. 
Uncorrelated atoms are well approximated by k = A^^/^ while highly squeezed states are obtained for k — j' 1. We 
consider the limit where A^ ^ 1 such that we can replace sums with integrals when calculating the moments of 
the angular momentum operators. This allows us to get analytical expressions for the moments of Jz in equation 



(22 1- (24 1. All the terms will decrease with growing k since AJz decreases for growing k, e.g., for a coherent spin state 



K = VN we have AJ^ '^ 0. For a fixed n all three terms will also decrease with growing n until a certain UmaxiK) 
is reached. At this point the uncertainty in J^ prevents us from information gained in further measurements. We 
believe that this effect is due to the finite probability of measuring a J^ with opposite sign than (J^), which spoils 
the measurement strategy. Note that rimax('*) grows with k since the width of J^ decreases with n. 



Of the three therm in Eqs. (22)- (24 1, the term in equation (24 1 is decreasing the slowest with n since this contains 
{1 — Jz/{Jz)) to the lowest power. In the following we therefore focus on this term. In order for the performance to 
be nearly Heisenberg limited we need the contribution from this term to be close to or smaller than then Heisenberg 
limit, and this put restrictions on the possible values of k. To determine the conditions for n we have numerically 
solved the equation ((1- J^/(J^))2"-2 j2/(j^)2^ ^ 1/7V2 gjj^^^g ^jj^i jg ^j^g Heisenberg limit of {6<^'i). This condition 
thus determines the parameters for which the noise from AJ^ is comparable to the Heisenberg limit. For TV in the 
range N = 10^ — 10^ we have identified the n for which the equation is fulfilled for the smallest k. The result is that in 
order for this term not to dominate the noise, n needs to grow with increasing TV to suppress the noise in J^, but the 
growth can be slower than k ^ log vN + 2 and requires n ~ 3 log A'' measurements. Thus if we choose k ~ log v iV + 2 



and n ~ 31ogA^ the noise terms in equation (22)- (24 1 will decrease as < l/N , which is the Heisenberg limit 



We now turn to the measurement noise, which consists of two parts. One is the accumulated back action contained 
((5J|„) (see equation (20)) while the other is due to the noise in the probe light (see the last term in equation 



(21)). For now we consider the accumulated back action. Using equation (15)-(17) we find that this is dominated by 



J27=i ^'i{^i){{^'t'tJx/{Jz)Y) and that the dominating terms from each measurement are 

z = l: \{sH5cl^,f){{l-lj{Jz)fjl/{Jz?)nl (25) 

z > 1 : \{{sH5cj,,) - 6M^){{l-Jz/{J.)r-^jVl/{Jz)')^1 (26) 

where we have assumed that the probe light has vacuum statistics such that (P) = {X) = {XP) = and (A^) = 
(P2) ^ 1/2. Again wc can get analytical expressions for ((l-Jj(J^))2jJ/(J^)2) and {{l-Jj{Jz)f^-'^JlJll{Jz)^) 
by using the Gaussian approximation. Motivated by the previous numerical calculations we set n = logyiV + 2, 
n ~ SlogA^ and by numerically evaluating the terms for A^ = 10"^ — > 10^ we find that for fli = iV^^^"+T all terms 
will be < l/A^^, i.e. at the Heisenberg limit. 

We now consider the part of the measurement noise that comes from the noise in the probe light. From equation 
(Il5])-([T7| and equation E^ we find that the dominating terms are YJlZi {{^- Jz/ {Jz)Y''~'^'){Pf) / {{Jz)'^^^)- Again 



by numerically evaluating the terms for A^ = 10"^ — > 10^ we get the scaling of the terms and we find that for 
K = \ogy/N + 2,nr^Z\ogN and ^^ = A"^+^ aU the terms will be < l/N"^. 

So far we have found that we can make the noise from the measurements and from AJ^ be < l/N"^, which is 



the Heisenberg limit. The limiting noise in ($„) is then the noise from AJy. From (20) we find that this has a 



contribution of {J^)/{Jzf = AJy/{Jzf. For the states \tp{K)) we find that AJy/{Jzf - k'^/N^. We thus get 

aj ~ (2/A^ + log ^/N/N)/^/jT for n = log y^ +2,n-^3logN and r^^ = A^"^+^ . Note that we are in the hmit of 
A^ ^ 1 and that a^ is in units of (7/(a;^r))^''^ (see article). For comparison the Heisenberg limit for the same Ramsey 
time is a^ = {l/N)/\/^T in the same units. Hence our results shows that near Heisenberg limited stability can be 
obtained with the adaptive protocol. 

We will also consider the performance of the adaptive protocol with uncorrelated atoms (k = ^/N) . In this case as 
for the conventional Ramsey protocol the stability will be limited by the projection noise and will be a-^ ^ N~^^'^ j \/^T . 
Furthermore we will compare with the optimized protocol of Ref.[5] where the stability is a-^ ^ A~^/^/V7T. 



NUMERICAL SIMULATION 

To verify the above findings we have simulated an atomic clock with a LO subject to both white and 1// noise. 
For both types of noise we have simulated the clock for atom numbers ranging from 100 to 10^. For N < 1000 we 
simulate the full quantum evolution during the measurements by bringing the input state \iP{k)) trough a Ramsey 
sequence and mixing it with a light state, which is assumed to be in vacuum. We pick the measurement outcome of 
P according to the corresponding probability distribution and subsequently update the state of the atoms for the 
next measurement etc. We denote this as 'full quantum simulation'. For A^ > 1000 we approximate the probability 
distributions of Jx,y,z with Gaussian distributions with moments calculated from \iP{k)) in the limit of A 3> 1 such 
that we can replace the sum over m with an integral. We denote this as 'Gaussian simulation'. 
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The adaptive measurement protocol and the feedback on the LO are simulated as described in equation (1) and 
below in our article. The feedback effectively locks the LO to the atoms, and lowers the noise level of the LO as 
shown in Fig. [4] where the noise spectrum S{f) of the LO is plotted against the frequency /. The noise spectrum is 
defined as S{f)S{f + /') = {5uj{f)Suj{f')) where Suj{f) is the Fourier transform of the frequency fluctuations (5w(t) of 
the LO. For a free running LO with white noise we use S{f) = 7 while for 1/f noise we use S{f) = 7^// where 7 is a 
parameter characterizing the fluctuations of the LO. The figure show that for high frequencies the locked LO has the 
noise of the free running oscillator but for low frequencies the LO is locked to the atoms and is limited by the atomic 
noise. Furthermore Fig. |4] shows how squeezing improves the stability of the clock by lowering the noise level of the 
locked LO more than for uncorrelated atoms. While the conventional Ramsey scheme works ideally for k ^ 14, the 
adaptive protocol allows for k '-^ 3 at the atom number N = 1000 used in the figure, and thus leads to an improved 
stability. 
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FIG. 4. (Colour online) Noise spectrum of the slaved LO for (a) white noise and (b) l//-noise. The noise spectrum is defined 
as S{f)5{f + f') = {5uj[f)5u){f')) where 5uj{f) is the Fourier transform of the frequency fluctuations 5uj{t) of the LO. The plots 
show how the feedback efi'ectively locks the LO to the atoms and the noise in the LO becomes limited by the atomic noise for 
low frequencies. Furthermore it is seen how squeezing lowers the atomic noise and hence the noise of the locked LO. The plots 
where made for A'^ = 1000 and "jT = 0.1 where the optimal performance of the conventional protocol is reached for a squeezing 
of K ~ 14. The adaptive protocol allows for k ~ 3 which further reduces the noise of the slaved LO. Note that the optimal 
performance of the adaptive protocol is reached at higher "fT. 



For a fixed Ramsey time the feedback strength (a) determines how long time the clock has to run before the LO 
is locked to the atoms. Since we simulate a clock running for a long but finite time there will be some remaining 
information from the last measurement results, which have not been fully exploited by the feedback loop. In our 
simulations we therefore do an additional phase correction to the LO after the final measurement. In principle the 
influence of the last few measurements could also have been reduced by running the simulation for a longer time, but 
by doing the correction we reduce the required simulation time. To flnd the required phase correction and obtain an 
expression for the stability of the clock we study the phase of the locked LO. At time t^ = kT the phase of the LO is 



54>{tk) 






(27) 



where Sujq (t) is the frequency fluctuations of the unlocked LO and AtUi is the frequency corrections applied at time 
ti. Using that Aw^ = —a6(j)e{ti)/T where 54>e{ti) is the estimated phase of the LO at time ti we can write 



fc+i 



5<t){tk) = dMtk)~aY,5(j>e{U), (28) 

where S(j>Q{tk) = J^^ SLL!o{t)dt. The mean frequency offset of the LO after running for a period r = ZT (/ ^ 1) is 



Suj{t) 






^filial correct 



(29) 
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where (pdnai correct IS the phase correction that we apply after the final measurement. Combining equation (281 and 
p9| we find that 0fi„ai correct = J2\^i ((1 - aY'^-ScfjeiU) + Sj=i "(1 - ay~^S(t>eitj)) wiU give the ideal performance. 
For this 0finai correct the mean frequency offset becomes 



SQ{t) 



1 ' 



T 



(30) 



We use this expression to determine the stability of the clock, which is given by (J-yir) = {(SLd{T)/uj)'^)^/^. Note that 
while the sum is over different time intervals we cannot in general determine the stability by looking at an interval 
independently since the phases are correlated for finite a or for correlated noise in the free running LO e.g. 1// noise. 
In our analytical calculations, however, we assume white noise and a -t> so that we can ignore the correlations and 
consider each Ramsay sequence independently. 

For the model investigated here, the stability increases with the Ramsay time but T is limited by two types of 
errors. For experiments or simulations running with a fixed Ramsay time there will always be a finite probability 
that the feedback loop jumps to a state with a phase difference of 27r (so called fringe hops 6 ) or that a phase jump 
that is large enough to spoil the measurement strategy occurs. For the adaptive scheme this happens for phase jumps 
> n while it happens for phase jumps > tt/2 for the conventional protocol. The reason for this is that the adaptive 
protocol is able to distinguish whether the phase lies in the intervals [0; 7r/2] or [tt/2; tt] since they will lead to different 
response when we rotate the state during the feedback. On the contrary the conventional protocol only has a single 
projective measurement and cannot distinguish in which of the two intervals the phase lies. In our simulations we 
see the phase jumps as an abrupt break down as we increase the width of the distribution of the acquired phase, e.g. 
(7^ = 7T for white noise. This break down is clearly visible in Fig. ^ (reproduced as part (a) of Fig ^ where we also 
show the similar plot for 1// noise (b)). Below we will investigate this breakdown in more detail. 




FIG. 5. Stability as a function of the Ramsey time {'jT) for white noise (a) and 1// (b) noise in the LO. The plots were made 
with A'^ — 10^. D, V are the best non-linear protocol of Ref. [5] while o, A are the adaptive protocol discussed in this article. 
The adaptive protocol allows for 7T ~ 0.3 and 0.2 for white and 1// - noise while the conventional protocol of Ref. [5] only 
allows for 'yT ~ 0.1 for both white and 1// noise. V, A correspond to uncorrelated atoms while o,n are the ideal choices of 
squeezing in the adaptive and conventional protocols respectively. The probabilistic nature of the errors that limits yT is seen 
in (a) where the stability for the adaptive protocol jumps back and forth for jT > 0.3 before it is definitely diminished. 



For a simulation running for a time r = IT, we have / samples of the accumulated phase of the LO during the 
Ramsey time T. Assuming that these samples have an independent Gaussian probability distribution with zero mean 
(this will be the case for white noise and is approximately true for 1// noise when the LO is locked to the atoms), 
the probability of all of these phases to be less than a critical value (a) is 



P{<a) 



1 — erfc 



W2a 



(31) 



where the variance of the distribution a depends on 7T. For large / and as a function of a this probability will drop 
abruptly from '-^ 1 to ~ around a certain a = amax- Defining amax to be the position where P{< a) = 1/2 we find 
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that 



Solving this equation gives 



1/2 = ( 1 - erfc ( -= ) ) . (32) 



(J Ytl (IT """^ J 1 \ 'J'-' / 

^ln(2/7r) + 2 ln(0 - 2 ln(ln 2) - hi(hi(2/7r) + 2 \n{l) - 2 ln(ln 2)) 
where we have expanded to first order in z = (1 — 2^^/') ^ ln(2)/L It is seen that the breakdown (cTmax) has a weak 



(logarithmic) dependence on I. Solving equation (32) with the Ihs. being equal to 0.95 and 0.05 we find that P{< a) 
drops from 0.95 to 0.05 within a window of ~ 2cr„ja3,/ln(^). Hence for large I the errors will appear very abruptly in 
the simulations. 

Ideally we should include correction strategies for the errors due to large phase jumps in our simulations (e.g. 
running with different Ramsey times would correct for fringe hops), but for simplicity we ignore this. This means 
that our simulations have a weak dependence on the number of steps we simulate, but since it is only a logarithmic 
correction we do not expect this to change our results significantly. Instead we find the upper limit of 'jT from the 
simulations plotted in Fig.[2J where I ~ 10®. 

When minimizing the stability for the adaptive protocol we numerically minimize in both the degree of squeezing, 
the number of weak measurements and the strengths of the measurements. For white noise we can use the expression 
for a.y in the limit of a ^ 1. For 1// noise we include the feedback on the LO as described in equation (2) and below 
in the article. Note that < a < 1 in order for the feedback to stabilize the LO [25j. We have used a ~ 0.1 in our 
simulations and do not expect our results to change significantly for a different choice of a. We also minimize the 
stability in the degree of squeezing for the conventional protocol of Ref. [5] for comparison with the adaptive protocol. 



